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ABSTRACT 

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period 
transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 ^m 
bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf 
stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+° °;^^^° (3.6 ^im) and 0.094%+° °^^^° 
(4.5 /im). The HAT-P-4b values are 0.142%+[]°}f^^ (3.6 /im) and 0.122%+°-[J}2%(4_5 ^^^^^ Tj^g 
planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and 
low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT- 
P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, 
along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3ct upper limits 
on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 fim) and < 0.085% 
(4.5 /im). We discuss these results in the context of the Spitzer secondary eclipse measurements of 
GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high 
eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. 
We place 3a upper limits on the |e coswj quantity (where e is eccentricity and uj is the argument 
of periapsis) for HAT-P-3b (< 0.0081) and HAT-P-4b (< 0.0042), based on the secondary eclipse 
timings. 

Subject headings: stars: planetary systems — eclipses - techniques: photometric 



1. INTRODUCTION 

There are about 177 confirmed exoplanets with orbital 
periods of less than 10 days and masses greater than 
0.1 Jupiter mass, which are often referred to as "hot 
Jupiters". About 145 of them transit their host stars^^. 
Hot Jupiters have high equilibrium temperatures, often 
over 1000 K, due to their proximity to their parent stars. 
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These planets make ideal targets for the Spitzer Space 
Telescope, which can directly detect their thermal emis- 
sion via time series photometry in the infra-red during 
secondary eclipse. 

The drop of total light from the planet-star system 
as the planet moves b ehind the star was first mea- 
sured independently by iCharbonneau et all (|2005f ) and 
iDeming et all (|2005l) . By measuring the eclipse depth 
in several photometric bands, one can construct a very 
low resolution emergent spectrum of the observed planet 
in the infrared (Charbonneau et al. 2008; Grillmai r et al.l 
I2008D . Comparison of similar spectra with models has 
suggested that there are two types of hot Jupiter atmo- 
spheres based on the presence or ab sence of temperature 
inversions in their upper layers (e.g., Knutson ct al.'200 
Machalck ct al. 2009; Todorov ct al. 2010; Bccrcr ct al 
lonFfPeming et al.ll26lH: iTodorov et al.ll2012l) . 

The mechanisms that create such inversions are not 
well established, but it is generally believed that planets 
with an inverted atmosphere have an additional opac- 
ity source in their atmospheres' up per layers, where the 
pressure falls below - 0.01 bar ()Burrows et al.l l200l 
.Fortnev et al.l [2008.) . In the past, this opacity source 
has been suggested to be gas phase TiO ( Hubcn-^ T" at al.l 
200l iBurrows et all [20071 [20081: iFortnev et al.. ,200631 



20081) . However, this molecule may form grains on the 



night sides and in cold traps deep in the atmospheres 
on the day sides of some hot Jupiters with inversions 
such as HD 209458b, where the pressure-temperature 
profil es are expected to c ross the condensation curve of 
TiO (ISoiegel et al.l [20091) . In addition, TrES-3 is suffi- 
ciently heated to have gas phase Ti O in its at mosphere, 
yet has no temperature inversion (jFressin et al. 2010.) . 
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while XO-lb is too cool to maintain it, but appea rs to 
have an inverted atmosphere (|Machalek et al.|[2008f) . On 
the other hand, a m ore recent 3D mixing s tudy of hot 
Jupiter atmospheres (jParmentier et alll2013| ) has shown 
that TiO may stay aloft in HD 209458b's atmosphere 
due to strong vertical mixing, if it forms grains no larger 
than several microns. According to this investigation, 
TiO in the form of larger particles should be depleted by 
the day-night cold trap. Thus, the role of TiO in the for- 
mation of temperature inversions remains controversial. 

There are at least two alternative explanations for 
the presence or a bsence of temperature inversions. 
iZahnle et al.l (|2009D suggested that sulfur compounds 
may account for the stratospheric opacity in some hot 
Jupiter atmospheres. Another hypothesis is that the 
presence or absence of temperature inversions is corre- 
lated with the magnetic a ctivity and r elated UV flux 
from the host star (.Knutson et al.l[2010( ). According to 
this idea, the increased UV flux received by planets or- 
biting active stars destroys the molecule responsible for 
the formation of temperature inversions. 

We can test these hypotheses by building up a large 
sample of hot Jupiters with secondary eclipse measure- 
ments and search for correlations with other system pa- 
rameters, such as stellar activity or the Ti abundance 
in the stellar photosphere. Although the Spitzer Space 
Telescope exhausted the last of its cryogen in May 2009, 
it still remains the best facility for carrying out these ob- 
servations. The observatory has continued operating at 
a higher temperature of approximately 29 K, cooled by 
passive radiation. Only two photometric bands are still 
available, the 3.6 and 4.5 /xm chann els of the Infra-Red 
Array Camera fIRAC. iFazio et al.l [2004'). Fortunately, 
measurements in these two bands are often sufficient to 
constrain the presence or absence of tem perature inver- 
sion in the upper layers of hot Jupiters (jKnutson et al.l 
[20T0I) . 

The energy budget is another aspect of transiting hot 
Jupiters that can be studied via secondary eclipse pho- 
tometry. These planets are expected to tidally evolve 
within ^ 1 Gyr of their formation, and their rotation 
periods should become eq ual to their orbital perio ds, as- 
suming zero eccentricity (jCorreia fc Laskadl2010l) . The 
transfer of heat from the day side to the night side can 
have strong influence on the emergent wavelength de- 
pendent flux of the day side of the planet, which is mea- 
sured directly by estimating the secondary eclipse depth 
at a given wavelength. Relating observations to atmo- 
spheric models can place constraints on the heat redis- 
tribution efiiciency an d the Bond albedo of the planet 
(jCowan fc Agoll r2011[) . The authors of this study sug- 
gest that planets with very high equilibrium tempera- 
tures have a narrow range of Td/Tp, where Td is the 
effective dayside temperature and Tq is the equilibrium 
temperature at the substellar point. This ratio is a mea- 
sure o f the redistribu t ion eff iciency of the atmosphere, 
which iCowan fc Agoll ()2011f) find to be uniformly low 
for the hottest planets, while cooler planets appear to 
exhibit higher range of heat redistribution efficiencies 
and/or albedos. 

In this paper we present Warm Spitzer time series pho- 
tometry of three short perio d transiting exoplanets - 
HAT -P-3b (iTorres et al.ll200l. HAT- P-4b fK ovacs et ahl 
I2007j) and HAT-P-12b (jHartman eral...2009l The phys- 



ical properties of the systems are liste d in Table [TJ In 
a previous paper (|Todorov et al.ll20T2l ). we focused on 
Warm Spitzer secondary eclipse photometry of three 
planets (X0-4b, HAT-P-6b and HAT-P-8b) with masses 
comparable to Jupiter's, orbiting F dwarfs. The mag- 
netic activity of stars as indicated by the Ca II HfcK 
line strengths is uncalibrated for effective temp eratures 
over 6200 K or under 4200 K (jNoves et al.lll984f ). Since 
the host stars of all three of these targets were above 
or close to the 6200 K boundary, our results did not re- 
liably test the correlation between magnetic activity of 
th e host star and the a tmospheric inversions proposed 
bv lKnutson et al.l (|2010[ ). which relies on Ca II HfcK line 
measurements for assessing the magnetic activity of the 
stars. 

Therefore, in order to better cover the available param- 
eter space, for this analysis we have chosen planets that 
orbit host stars with effective temperatures within the 
4200 — 6200 K range. Our targets are less massive than 
Jupiter (the least massive, HAT-P-12b, has a mass simi- 
lar to Saturn's), but have radii similar to Jupiter's. In ad- 
dition to the secondary eclipse depth measurements, we 
combine their timings with the most precis e ephe merides 
available (HAT-P-3, HAT-P-4: iSada et al.l (poTl . HAT- 
P-12: this paper) to constrain the orbital eccentricity of 
the planets. 

For this analysis, we update the code developed by 
iTodorov et al.l (|2012[ ). to make use of full array mode ob- 
servations (in the cases of HA T-P-4 and IIAT-P-12) and 
variable aperture photometry (jLewis et al...2013 ) . We es- 
timated the uncertainty of the eclipse depth and timing 
with the Markov Chain Monte Carlo and prayer-bead al- 
gorithms. It is in principle possible to use this code to 
analyze all available Spitzer secondary eclipse data sets 
in a single investigation, but we find that every data set 
has peculiarities that need to be addressed on individ- 
ual basis. We analyze and report eclipse results for hot 
Jupiters in groups of three because this allows for the ef- 
ficiency of analyzing data for multiple exoplanets, while 
retaining the ability to cope with the quirks of individual 
data sets. 

In Section [2] we present the Spitzer observations and 
our photometry procedures. The details of the time se- 
ries analysis and the uncertainty estimation are presented 
in Section[3l We discuss our results in the context of pre- 
vious studies in Section |4l 

2. OBSERVATIONS AND PHOTOMETRY 

2.1. Secondary Eclipse Observations with Spitzer 

We observed the secondary eclipses using the Infra- 
Red Array Camera (IRAC) on the Spitzer Space Tele- 
scope in both the 3.6 /xm and the 4.5 /^m channels. The 
observations on the HAT-P-3 system in both pass-bands 
utilized the IRAC subarray mode, resulting in 32 x 32 
pixel (39"x39") images, centered on the planet's host 
star. The HAT-P-4 and HAT-P-12 data are in full ar- 
ray mode, where each image is taken has the full IRAC 
resolution of 256 x 256 pixels (5.2' x 5.2'). 

The subarray mode images are stacked in FITS data 
cubes that contain 64 exposures taken in a sequence. 
Our HAT-P-3 observations have effective exposure times 
per image of 1.92 s in both pass-bands. The full array 
data have effective exposure times of 4.4 s (HAT-P-4) and 
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10.4 s (HAT-P-12) per image respectively in both wave- 
lengths. Each secondary eclipse observation of HAT-P- 
3b lasted for 7h 42min resulting in 13,760 images (215 
data cubes). The observations of HAT-P-4b covered 7h 
38min, resulting in 3,871 images in each wavelength, 
while those of HAT-P-12b lasted almost as long - 7h 
37 min - but resulted in only 2,097 images per pass-band, 
due to the longer exposure time. 

Complete information about the time span of the ob- 
servations is presented in Table O 

2.2. Photometry and Time Information Extraction 

Our time-sta mp extraction r o utine is very similar to 
the one used bv iTodorov et al.l ()2012D . In all data sets, 
we perform the extraction on the Basic Calibrated Data 
(BCD) files produced by version SIS. 18.0 of the Spitzer 
pipeline. For timing of the photometric points, we rely 
on the MJD.OBS keyword in the FITS headers, cor- 
rected to indicate the mid-exposure time. We convert 
this time stamp, given in modified Julian date, based 
on the Coordinated Universal Time (UTC) standard 
(MJDuTc), to Barycentric Julian Date based on the Ter- 
restrial Time standard (BJD tt) using Jason East man's 
IDL routine get_spitzer_bjd (jEastman et al.l [20101) . We 
prefer to use the TT standard rather than UTC, be- 
cause the latter is discontinuous and has leap seconds 
introduced occasionally. At the time of our observations, 
TT « UTC+66.184S. 

We convert the pixel intensities from MJy/sr to elec- 
tron counts using the the information provided in the im- 
age headers, in order to be able to estimate the Poisson 
noise of the photometry. The data are filtered for ener- 
getic particle hits by following each pixel through time. 
This is done in two passes, first flagging all pixels 8 a 
or more away from a boxcar median through time with 
width 5. Their values are replaced with the local boxcar 
median value. In the second pass, all values more than 
4(7 away from the boxcar median through time (again 
with width 5) are flagged and their values are replaced 
with the local median value. For the HAT-P-3 data, this 
procedure is performed separately for each data cube in- 
stead of for the whole time series in order to avoid pixel 
rejection due to the sharp changes in background that 
occur between the last frame from one cube and the 
first frame f rom the next. This effect is similar to the 
one seen bv iDeming et al.l (|2011l ). The fraction of cor- 
rected pixels is about 0.53% (HAT-P-3, at 3.6 /xm), 0.12% 
(HAT-P-3, at 4.5 Aim), 0.13% (HAT-P-4, at both chan- 
nels), 0.22% (HAT-P-12, at 3.6 /xm) and 0.24% (HAT-P- 
12, at 4.5 /^m). 

We estimate the background flux from the full-array 
images by creating a histogram of all pixel values for each 
frame and fitting a Gaussian function. We correct for 
contamination by field stars and energetic particle hits 
by fitting only to the central portion of the histogram 
and excluding the regions that account for the high pho- 
ton count pixels. We test for variable background across 
the full-array images by measuring the background based 
on 30x30, 50x50, and 100x100 pixel boxes centered on 
the star-planet system (which was always positioned in 
the center of the array). Removing the background from 
the photometry in this manner produces at most 0.7% 
change in the photometric scatter around a running me- 
dian of width 20 of the photometry compared to data 



reduction using the whole array to estimate background, 
in all data sets. This is marginal, therefore, we elect to 
maximize our background determination precision and 
use the whole array to determine the background lev- 
els. Since the subarray mode observations of HAT-P-3 
result in only 32 x 32 pixel images, we exclude a 5 x 5 
pixel square centered on the star from the histograms 
for these data, in order to avoid biasing the background 
estimation to higher values. 

In order to locate the centroid of the stellar point re- 
sponse function (PRE) we experiment with fitting a two- 
dimens ional Gaussian f unction to the core of the stellar 
image (lAgol et al.ll2010D and with flux- weighted cen troid- 
ing (e.g.. lCharbonneau et alll2008l: iKnutson et al.| [2008'). 
We find that in most cases there is an improvement in 
the standard deviation of the residuals that remain after 
we subtract our best fit model from the photometry, if 
we adopt flux-weighted centroiding. The only exception 
is the HAT-P-3b 4.5/im data set, where the difference in 
the resulting scatter is marginal (~ 0.1%). Therefore, we 
adopt flux-weighted centroiding for all data sets in this 
analysis. 

We experiment with two different photometry ap- 
proaches. First, we perform aperture photometry, using 
the IDL routine aper^'^ , varying the aperture radius in 
increments of 0.5 pixels between 1.5 and 6.5 pixels. We 
select the best photometry aperture radius by measuring 
the true scatter it produces around a boxcar median of 
the raw light curve with width 20. The scatter is not 
strongly correlated with aperture radius. For the HAT- 
P-3 subarray mode data, the minimum scatter is found 
at 2.5px (3.6 /im) and 3.0px (4.5 /im) pixel radii. For the 
full array mode data we find that the smallest scatter 
occurs at pixel radii of 4.0px (HAT-P-4, 3.6 /im), 3.0px 
(HAT-P-4, 4.5 /^m), 5.0px (HAT-P-12, 3.6 /^m) and 2.5px 
(HAT-P-12, 4.5 Aim). 

Separately, we perform photometry on all data by us- 
ing time- variable aperture. For each image, we estimate 
the noise-pixel parameter ([Mighelll 120051 : iKnutson et al.l 
2012; Lewis et al. 2013), which is a measure of the width 
of the stellar point spread function (PSE). It is defined in 
Section 2.2.2 of the 5piteer//i?^C instrument handbook 
as: 

/3 = §^, (1) 

where li is the intensity detected by the i^^ pixel. The 
noise pixel parameter, /3, is proportion al to the ful l- 
width-half- maximum of the stellar PSE (iMighelll 120051) . 
For each image, we calculate the photometric aperture 
radius: 



(2) 



where b is a scaling factor and c is a constant. In each 
image frame, we measure the flux used to determine /3 
using circular aperture radii between 1.0 and 6.5 pixels. 
If any part of a pixel falls within the aperture radius, it 
is fully included in the /? calculation. For each of these 
values, we vary h and c in steps of 0.05. We fit each 
of the resulting light curves with our "systematics-and- 
eclipse" model and find the combination of photometric 

^"^ http:/ /idlastro. gsfc.nasa.gov/ 
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parameters that yields the smallest standard deviation 
of the residuals. 

For the 3.6/xm data sets, the variable aperture photom- 
etry approach yields lower residual scatter values (be- 
tween 3 and 7%), reduces the amplitude of the periodic 
flux oscillation in the raw data by about 50% and re- 
duces the levels of correlated noise in the residuals after 
the best fits to the data are subtracted from the photome- 
try. This is particularly evident in the HAT-P-12b 3.6/im 
light curve, where fixed aperture photometry yields sig- 
nificantly different eclipse depths (between —0.07% and 
0.07%, which is above our final 3a eclipse depth limit) 
based on the radius of the photometric aperture and the 
method for decorrelation of the intra- pixel effect. 

In the 4.5/xm light curves, the residual scatter values 
are consistently, but marginally smaller for the variable 
aperture than for fixed aperture (between 0.1 and 1%). 
There is also a slight improvement in the residual corre- 
lated noise levels in the light curves after decorrelation. 
Therefore, we adopt variable aperture in the final analy- 
sis for all data sets. The photometry parameters we used 
are summarized in Table [2] The raw photometry during 
secondary eclipse for the three planets is presented in 
Figure [TJ 

3. DATA ANALYSIS 

3.1. Ephemerides 

The data analysis routines fit the eclipse models to 
the light curves as a function of orbital phase. There- 
fore, the best available estimates of the ephemerides of 
the planets are needed in order to calculate the orbital 
phase with minimum uncertainty. For HA T-P-3b and 
HAT-P-4b, we use the ephemerides given in iSada et all 
([2012). The ephemeris provided in that paper for HAT- 
P-12b, how ever, does not in clude the transits observed 
recently by iLee et ahl (|2012| ) , and vice versa. We com- 
bine the transit timings provided by both groups with 
other transit timings from the literature (Table |4]). The 
resuhs are Tq = 2454187.85559 ± 0.00011 in BJDtt and 
P = 3.21305929 ± 0.00000034days. The offsets in min- 
utes between the best fit ephemeris and the transit tim- 
ings are shown in Figure [21 and they are consistent with 
an unperturbed orbit. 

3.2. Secondary Eclipse Fits 
3.2.1. Data Examination 

The secondary eclipse depth is only measurable in 
Spitzer time series photometry after careful removal of 
any instrumental effects. We normalize the light curve 
so that the mean brightness of the target system during 
eclipse is unity, corresponding to the light of the star only. 
We remove any data that have high backgrounds or are 
outliers (9 frames for HAT-P-3 at 3.6/im, 5 frames for 
HAT-P-3 at 4.5 /xm, 22 frames for HAT-P-4 at 3.6 ^m, 
42 frames for HAT-P-4 at 4.5 fim, 75 frames for HAT- 
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A^-ol et al. 2010; D eming et al.ll2011[ 

Todorov et al.u20l3) . we find that the 57*'' frame in each 



Previous time series photometry with the 
Spitzer/IRAC instrument, both during the cryo- 
genic and warm missions, has revealed two different 
transient instrumental effects. First, it often takes tens 
of minutes for the target star's positi on on the detector 
to stabilize (e.g., [Anderson et al. 20Tl[ ). This initial 
position instability causes apparent changes in intensity 
because it involves portions of the detector that have 
different intrapixel sensitivity variations than the ones 
used the most during the observa tion . Th e second 
effect, seen by e.g., Canipo et al. (2 0111 ): iDem ing et aL| 
(|2011l ): [ibdorov et all (|2012l ). causes the apparent 
brightness at the start of some observations at 3.6 /im 
to increase or decrease in an exponential-like manner 
before stabilizing, without correlation with the position 
of the stellar image on the detector. This behavior is 
similar to that observed in the longer-wavelength IRAC 
arrays and is believed to be due to char ge-trapping (e.g., 
iKnutson et all l2007t lAgol et al.l I2O1O0 . The simplest 
way of correcting for these effects is to clip the initial 
portion of the time series. Therefore, we discard the 
initial 48min 42 sec from the HAT-P-3 data at 4.5 /im, 
corresponding to 1423 frames taken before orbital phase 
of 0.44. We do not find it necessary to clip the HAT-P-3 
at 3.6 //m data, or the HAT-P-4 and HAT-P-12 time 
series in either wavelength. The corrected light curves 
are shown in Figure [3] 

3.2.2. Initial Fitting Procedure 

Our procedure for determining the most suitable 
systematics and ech pse model is similar to that of 
iTodorov etaP (poTl ) - we assume a central phase of the 
eclipse and perform a simultaneous linear regression fit 
for all free parameters to the unbinned light curve. We 
make incremental increases to the initial assumed phase 
to scan a wide range of possible central phases of the 
eclipse, and fit a new model to the data, after each step 
in phase. We track the values of the best regression 
fits as we make the scan. The step size in phase is 10"^ 
for all data sets, and we cover the phase intervals be- 
tween 0.48 and 0.52. We experiment with larger phase 
ranges, especially in the case of HAT-P-12b, but without 
an improvement in the results. 

We employ a computational model describing the sys- 
tematic and astrophysical effects observed in our data 
sets similar to the one used by ITodorov et a l. (2012 



Similarly t o previous studies (^e.g.. [ Knutson et al.ii200 



data cube in the subarray mode exhibits a relatively high 
background value, and we exclude all these images from 
the analysis (215 frames in each HAT-P-3b band). 



i Beerer e t al. 2011; Deming et al.il2oil[ ). we find a correla- 
tion between the X and Y position of the stellar image on 
the pixels and t he me asured brightness of the star. As in 
ITodorov et all (poTl ). for a given eclipse central phase, 
we adopt a model of the systematics with a quadratic 
dependence between X and Y in intrapixel coordinates 
without cross-terms and the intensity and a linear ramp 
with time: 

I{t) ^at + biX + b2X^ + ciY + C2Y^ + Io + eiM, (3) 

where I(t) is the intensity as a function of time, t is time 
in units of phase, X and Y are the positions of the stel- 
lar centroid on the pixel in the x and y directions, M 
eclipse shape model, and a, 6i, 62, ci, C2, /o (the ordinate 
axis intercept) and ei (the eclipse depth) are the free pa- 
rameters. We also experimented with adding the noise 



Warm Spitzer Photometry of Transiting Hot Jupiters HAT-P-4b, HAT-P-4b and HAT-P-12b 



5 



pixel parameter, /3, as a third dimension in the spatial fit, 
but we noticed no improvement in the quality of the fits. 
Thus, we have not included /3 in the fit. Multiple previ- 
ous studies hav e also settled on quad ratic X and Y decor- 
relati o n fe.g.. [Charbonneau et aLI 1^ 008; Kn utson et al.l 
200a iChristiansen et al.l l2010t lAj^derson et al.l 12011 



Cochran et al.l 120111: iDemorv et al.l 120111: iDesert et alj 
2011f) . We experiment by setting all combinations of a, &2 
and C2, to and examine the residuals, but we find that 
the lowest residual scatter occurs when all parameters 
are left free (as expected). 

We, therefore, attempt to use the Bayesian information 
criterion test (BIC) to determine the optimal number of 
free parameters. We find that for HAT-P-3b (both wave- 
lengths) and HAT-P-12b at 3.6 /^m setting &2 = results 
in minimum BIC. The BIC for HAT-P-4b at 4.5 ^m is 
minimized by setting all parameters free. The minimum 
BICs for HAT-P-4b at 3.6 /xm and HAT-P-12b at 4.5 ^iuy 
are found when &2 and C2, are set to 0. However, for these 
data sets, this results in red noise with large amplitude in 
the residuals, and an unrealistically large eclipse depth 
value for HAT-P-4b at 3.6 /im. Setting 62 and C2, free 
minimizes the red noise amplitudes in these light curves. 
For all other data sets, regardless of the choice of pa- 
rameters, the eclipse depth and central phase values are 
within one sigma of each other. We experiment with 
adding higher order terms to Equation [3l but this does 
not lead to improvements in the red noise reduction or 
the BIC values, and to only to marginal improvements 
in the scatter of the residuals. Therefore, we conclude 
that the BIC test is not ideally suited for the correction 
of the systematic noise in our data. 

To further motivate this, suppose we found two curves 
that accounted for the intra-pixel effect, and their best- 
fits were essentially the same curve (i.e. they lay on 
top of each other), but one used many more parame- 
ters. We certainly would not be justified in identifying 
the additional parameters with the physical properties of 
the detector (lacking a priori information on the physics 
of the effect). However, as far as removing the effect 
from the photometry, we could use either curve because 
they would produce the same decorrelated photometry, 
although their BIC values might be very different. Thus, 
the BIC is not necessarily relevant, because we're not 
seeking information about the detector. The BIC would 
be valuable in distinguishing various decorrelation mod- 
els if they did not overlap each other, and one curve pro- 
duced a significantly smaller value. However, in our 
data adding cubic and higher terms to the decorrelation 
polynomial produces only marginal decreases in the x^. 

Thus, we elect not to use the BIC to determine the 
optimal number of free parameters and choose to keep 
all parameters in Equation [3] free, while not adding any 
higher order terms. 

The measured eclipse depths for the HAT-P-12 data 
sets are consistent with zero. We experiment by remov- 
ing various portions of the data at the start and the 
end of the observations. The resulting fits have cen- 
tral phases covering the whole explored range, and the 
eclipse depths take small negative or positive values. We 
conclude that our photometric precision is insufficient to 
detect the eclipse in these data sets and we place upper 
limits on its depth. 



For completeness, we experiment with substituting 
the quadratic dependencies of intensity on X and Y 
with a weighting function, as described bv lBallard et al.l 
()2010a| ). which essentially multiplies each photometric 
point by a "weight" , dependent on the X and Y positions 
of the stellar image on the detector. The apparent bright- 
ness of the star is smoothed as a function of X and Y on 
the IRAC array, with smoothing widths Cx and Cy. The 
in-eclipse data are excluded to prevent the decorrelation 
from removing it as a systematic effect. We optimize the 
smoothing widths to minimize the scatter in the residu- 
als of the data after subtracting the best fit model. For 
the 3.6 /im data we settle on CTx = 0.00425,0.00638, and 
0.00723 and ay = 0.00723, 0.00638, and 0.01064 for HAT- 
P-3b, HAT-P-4b and HAT-P-12b, respectively. For the 
4.5 /im data, we use cr^ = 0.00894,0.03404, and 0.03404 
and CTy = 0.00936, 0.02085, and 0.03574, again for HAT- 
P-3b, HAT-P-4b and HAT-P-12b, respectively. In this 
case, the BIC is not an adequate measure to compare the 
resulting fits, since this is a non-parametric correction of 
systematic effects (for more details on the applicability of 
the BIC, see e.g., iStevenso n et al.l[2012t ). Therefore, we 
compare the mean standard deviations of the residuals 
after subtracting the models from the photometry. We 
find that the scatters produced by weighting function fits 
are comparable or slightly larger than the ones from our 
polynomial fits, and using the weighting function decor- 
relation does not improve the correlated noise removal. 
Therefore, we elect to use quadratic fits in X and Y for 
all data. 

3.2.3. Best Parameter Values and Uncertainty Estimates 

In order to determine the best parameter values and 
their uncertainties, we utilize two approaches - Markov 
Chain Monte Carlo (MCMC) and prayer-bead Monte 
Carlo (PBMC). We first discuss our implementation of 
these algorithms, and then focus on our best fit param- 
eter determination approach. 

We implement a Markov Chain Monte Carlo (MCMC) 
code in order to quantify the uncertainties on the eclipse 
depths and central phases tha t we measu re. We follow 
the recipe suggested by ,Ford| ()2005l I2006D , taking steps 
in one parameter at a time, and drawing the steps from a 
Gaussian probability distribution. We perform 10^ iter- 
ations per free parameter. Here, we add central phase as 
a formal free parameter. Before running the main chain, 
we run several shorter chains to optimize the most likely 
step size for each parameter result in acceptance rates 
between 35% and 55%. A typical acceptance rate for 
all parameters and all data sets is ^ 45% , which is near 
the ideal rate suggested by iFordI (|2006L and references 
therein). We run the MCMCs for 7 x 10'^ steps. 

Using the steps from the Markov chains, we create his- 
tograms of the eclipse depth and central phase values 
during a given run. These have shapes close to Gaus- 
sian (for the eclipse depth histograms, see Figure 01) and 
we estimate the uncertainties on the astrophysical pa- 
rameters by calculating their standard deviations from 
the best fit values. MCMC histograms tend to be close 
to Gaussians even for data sets dominated by correlated 
noise, since this algorithm assumes white noise. 

The MCMC fails to converge on a solution for the two 
HAT-P-12 data sets, due to the undetectable eclipses. 
Hence, we report the upper limit for the eclipse depths 
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based on the regression analysis uncertainty, assuming 
central phase of 0.50007. This number includes the light 
travel time delay assuming ecos(w) = 0, but does not in- 
clude any apparent delay that may be due to the hottest 
point on the planet trailing behind the substellar point 
along the orbit. This ef fect causes a typically small de- 
lay, a bout 20 -30s (e.g., iKnutson et"aLll2007l: ' lAgol et al.l 
[2010h . 

Siin il arly to, e.g.. iDesert et al.l (|2011[ ); iDeming et al.l 
(|2011| ): iTodorov et al.l (|2012| ). ^ obtain an estimate of 
the systematic uncertainties due to correla ted red noise 
by pe rforming a "prayer-bead" analysis (|Gillon et al.l 
l2007a[) . In this method the residuals of the best regres- 
sion fit are shifted right by one frame (last frame becomes 
first) and added back to the best fit model, thus creat- 
ing a simulated data set with the red noise preserved. 
We perform a fit to the simulated data set using the 
algorithm described in Section 13.2.21 We record the re- 
sulting eclipse depth and central phase, and simulate an- 
other data set by shifting the residuals of the original 
best fit again. The histograms of the resulting eclipse 
depth distributions are presented in Figure |4l They are 
non-Gaussian, as expected for red noise dominated data 
sets, but begin to approach a Gaussian for data where 
the white noise dominates. 

The smallest value of a fit represents the best fit 
value only in data sets dominated by Gaussian noise. 
The Spitzer light curves, however, are often dominated 
by red noise, which is not correlated in time with the 
astrophysical signal. Therefore, for a red noise domi- 
nated light curve, any of the data sets simulated in the 
PBMC simulation run could have been the observed data 
set. Hence, for both eclipse depth and central phase, we 
elect to report our best fit value to be the median in the 
histogram of a parameter from the MCMC or PBMC 
fit simulations, whichever yields the larger uncertainty 
range. We adopt the corresponding uncertainties (Fig- 
ureU). If the data sets are dominated by Gaussian noise, 
our approach reduces to adopting the parameter values 
that yield the smallest x^- 

The best fit eclipse depths from the original data are 
very close to the medians from the Monte Carlo runs. 
For the IIAT-P-3 data, the original data best fits are 
0.108% and 0.096% versus Monte Carlo median eclipses 
of 0.1121° ;;;^^ and Omitoml at 3.6 and 4.5 /im, re- 
spectively. In both cases, this is a difference of about 
0.13cr. The HAT-P-4b original data eclipses are 0.142% 
and 0.124% versus 0.1421° °^^ and 0.122t°;°}^ from the 
Monte Carlo runs at 3.6 and 4.5 fim, respectively. The 
difference between these values is Oct (3.6 /xm) and 0.17ct 
(4.5 /im). Thus, the exact choice of "best values" has no 
impact on our results. 

We summarize the Monte Carlo results we adopt as 
final in Table [5] The uncertainties in central phase in 
that table include an uncertainty contribution from the 
ephemeris, but this is small (~ 1% of the total timing 
uncertainty). The eclipse central times in BJDtt are 
independent of any uncertainty in the ephemeris. 

4. DISCUSSION 

4.1. Comparison to Models 

We compare our eclipse depth measurements to two 
sets of models, by [Burrows et al., (|2007> i2008< ) and 



iFortnev et all (|2005l l2006allbl . l2008l) . in order to bet- 
ter understand their implications for the thermal struc- 
tures and heat transport efficiency of the planetary at- 
mospheres (Figure [5]). The chemical equilibrium and 
opacities i n the B urrows models are based on studies by 
,Burrows fc SharpI fl999. ) and Sharp fc Burrows {lOOlj, 
respectively. The parameters in these models that are 
relevant to us are Kabs, the absorption coefficient of the 
unknown stratospheric absorber, and Pn, the heat re- 
distribution parameter, which describes the amount of 
stellar flux transported from the day side of the planet 
to its night side. Kabs has units of cm^g~^, while Pn is 
unitless and varies between (no redistribution) and 0.5 
(complete redistribution). 

In the Fortney models, essentially, only the heat re- 
distribution efficiency, /, is a free parameter. The high 
altitude TiO and VO at equilibrium abundances serve as 
absorbers causing temperature inversions when needed 
to explain the data. The heat redistribution efficiency 
varies between / = 0.25 (indicating that the flux is evenly 
distributed over the whole planet) and / = 0.67 (corre- 
sponding to no flux redistribution at all, even within day- 
side regions of different temperature). The value / — 0.5 
signifies that fiux is evenly redistributed on the day-side 
of the planet, but no heat leaks to the night side. The 
atmospheres in the Fortney models have solar compo- 
sitions, except the / = 0.6 IIAT-P-12b model (shown 
in red in Figure [5]), which has a 30x solar metallicity. 
Neither the Fortney nor Burrows models account for the 
presence of clouds or disequilibrium chemical processes. 

"Fitting" these models in the mathematical sense is not 
practical. In the Fortney models one can vary / and the 
presence or absence of TiO and VO, while the Burrows 
models have Kabs and Pn as free parameters, but neither 
set of models is intended to include an algorithm to ad- 
just these parameters according to the data as part of the 
model calculation. The best approach available with the 
current state of the model codes is to compute a range of 
models for a given planet covering the parameter space, 
and then to select manually the one that accounts best 
for the data. For HAT-P-3b, we visually examine Bur- 
rows models with Kabs between and 0.1 crri^ g^^, and 
Pn between 0.1 and 0.3: for HAT-P-4b we look at ftabs 
between and 0.2 cm^ g"^ , and Pn between 0.1 and 0.3; 
and for HAT-P-12b we experiment with Kabs between 
and 0.1 cm? g^^ , keeping P„ at a moderate value of 0.3. 
We also examine Fortney models with and without TiO 
and VO in the upper layers of the atmosphere, with / be- 
tween 0.25 and 0.6. Below, we discuss models that match 
the data well, and our conclusions for the atmosphere of 
each planet. 

4.2. HAT-P-3b and HAT-P-^h 

HAT-P-3 displays a slightly enhanced level of chro- 
mosphcric activity, with calcium HfcK a ctivity index 
log(R^K) = -4.904 (jKnutson et al.ll20Tol ). This could 
explain some of the red noise evident in the light curves 
and PBMC runs for this planet (Figures |3| and |4|) . 

For HAT-P-3b (top panels in Figure O, we adopt the 
Burrows model with Kabs = 0.1 cm^ and Pn = 0.1, 
indicating an atmosphere with a temperature inversion 
and modest heat redistribution. We find that the atmo- 
sphere of IIAT-P-3b is matched by a Fortney model with 
/ = 0.6. This planet, like HAT-P-12b, is cool enough 
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that the presence of TiO and VO has practically no ef- 
fect on the shape of the spectra, and the Fortney mod- 
els cannot be used to distinguish between inverted and 
non-inverted temperature profiles. HAT-P-3b's eclipse 
depths are matched equally well by the almost identical 
inverted and non-inverted spectral models. 

Both sets of models, however, agree on low redistribu- 
tion ef ficiency, in apparen t contradiction of the hypoth- 
esis of iPerna et al.l ()2012[ ) that hot Jupiters with irradi- 
ation temperature, Tirr ^ 2000 K have efficient flux re- 
distribution. In their convention, Ti^ — Toff (R*/a)^/^, 
where Teg is the effective temperature of the star, 
is the radius of the star and a is the semimajor axis. 
For HAT-P-3b, Ti„. = 1600 K. Our result does not con- 
tradict (but do e s not support either) the hypothesis by 
iCowan fc Agoll ()2011| ). who suggest that a wide range 
of redistribution efficiencies are possible for planets with 
Te=o < 2400 K, since T,=o = 1300 K for HAT-P-3b. 

HAT-P-4 is chromospheric ally quiet, with an ac tivity 
index of log(R^K) = -5.082 (jKnutson et al.ll2010t ). and 
we find no significant perturbations to the light curves 
other than the eclipses. This planet is hotter than HAT- 
P-3b and HAT-P-12b and therefore has deeper eclipses, 
despite the fact that its host star is more luminous than 
the host stars of the other two. 

The HAT-P-4b models are presented in the middle 
panels in Figure [H We find that a Burrows model 
that describes the data well has Kabs — 0.2 cm^ g^^ and 
Pn = 0.1, corresponding to an inverted atmosphere with 
inefficient heat redistribution. The Fortney models with 
and without TiO and VO absorption in the upper atmo- 
sphere and / = 0.5 also seem to be close to the obser- 
vations. As in the HAT-P-3b case, the models appear 
to be ambiguous about any temperature inversions, but 
agree that the planet's atmosphere has moderate to low 
efficiency in redistributing heat to the night side. 

For HAT-P-4b, Ti„. = 2400 K. Therefore, the ineffi- 
cient flux redistribution we observe is consistent with 
the idea that planets with Tii,. above 2200-2 400 K should 
have little flux transfer to their night sides (jPerna et al.l 
[2012; Cowan & Agol 20111). 

iKnutson et al.. (,2010i ) hypothesize that while chromo- 
spherically active stars have planets with non-inverted 
atmospheres, planets around quiet stars have tempera- 
ture inversions. The difference between the activity in- 
dices of the host stars HAT-P-3 and HAT-P-4 is mini- 
mal (Alog(RHj^) = 0.1), and both fall in the tentative 
border region where hot Jupiters can have inverted or 
non-inverted atmospheres and so predictions for the pres- 
ence or absence of an inversion are difficult. Due to this 
and the ambiguity of the models of the planetary atmo- 
spheres, we cannot make claims that support or contra- 
dict this idea based on our data on these two planets. 

4.3. HAT-P-12h 

In principle, it is possible that we have failed to de- 
tect the HAT-P-12b eclipses due to a rel atively large or- 
bital e ccentricity. The discovery paper bv lHartman et aP 
(|2009( ) fixes the eccentricity at 0, which is what we 
have assumed in our determination of the upper limit 
of the eclipse depth. However, their initial fit results 
in a best value of e cos a; = 0.052 ± 0.025, which is in- 
significant. We estimate that we would have detected 



eclipses centered between phases of 0.45 and 0.55 (the 
range of our data in units of orbital phase) , correspond- 
ing to |e cosa;| < 0.08. This is only a bout la from 
the insignificant value bv iHartman et ah] (2009). Using 
the Exoplanet Or bit Database (http://exoplanets.org, 
iWright et all 1201 Ih . we find that out of the 188 known 
transiting planets with periods less than 10 days, only 
23 {'^ 12%) have orbital eccentricity over 0.08. Thus, 
we conclude that it is possible that we have missed the 
eclipse due to eccentricity larger than ~ 0.08, but that 
weak eclipses are a more likely explanation. 

Despite the formal non-detection, the HAT-P-12b 
4.5 /im light curve in Figure [3] appears to the eye to 
contain an eclipse. Allowing for a quadratic out-of- 
eclipse variation (due to, e.g., stellar spots or phase- 
variation), then we find a non-zero eclipse depth of 
- 0.049% ± 0.021% (a - 2.5cr resuh; still not a de- 
tection). This result is well within one sigma of the 
best value in Figure SI but is in better agreement with 
the non-inverted Burrows model. However, none of the 
other light curves requires anything but a flat out-of- 
eclipse baseline. Using a quadratic curve to allow for 
phase variation will naturally cause any eclipse in the 
data to look deeper, since the assumed maximum of the 
quadratic curve is near mid-eclipse, above what is ex- 
pected for a flat baseline. On the other hand, the ap- 
pearance of an eclipse and a phase curve variation is 
likely to be caused by instrumental red noise, given the 
lack of similar signatures in the other two 4.5 /xm curves. 
HAT-P-3b and HAT-P-4b are hotter than HAT-P-12b 
and have deeper eclipses, so any phase curve variability 
should be higher for them. In addition, the HAT-P-12 
host star ha s a Ca II HfcK activ ity index, log(R^j^), of 
about —5.1 (jKnutson et al.l[20Tol ). so stellar spots should 
be a less important factor than in the ligh t curve than 
they are, e.g., in the Sun (log(R{jj^) = — 4.9. lNoves et al.l 
[1981 . Hence, we conclude that a quadratic out-of-eclipse 
variation does not improve our results. 

We are only able to place upper limits on the secondary 
eclipses of HAT-P-12b, (assuming |e cosa;| = 0), how- 
ever, we can still compare these results to atmospheric 
models (lower panels in Figure[5|) . Both the inverted Bur- 
rows (O.lcm^g"^, Pn = 0.3) and the Fortney (/ = 0.5 
and 0.6, with 30 times solar metallicity for the latter, 
with the presence or absence of TiO irrelevant to the 
shape of the spectrum, due to the low temperature of 
the planet) models are poor flts to the data. The non- 
inverted Burrows model (/Cabs = 0, Pn = 0.3) appears 
close to the Sa eclipse depth upper limits, consistent with 
a lack of temperature inversion. 

The non-detections of thermal emission at both Spitzer 
channels, may suggest that heat is relatively efficiently 
redistributed to the night side of HAT-P-12b and/or 
that the albedo is high. However, it is difficult to claim 
this with any certainty, given the poor match that the 
models provide for the HAT-P-12b measurements. For 
this planet, Tjrr = 1350 K, but its surface gravity is ^ 
5.7ms^^, significantly less than the value of 10 ni s ~^, as - 
sumed for "typical" hot Jupiters bv IPerna et al.l (|2012D . 
For these two reasons, we do not consider the possible 
high heat redistributio n efficiency o f HAT -P-12b to be 
evidence in favor of the IPerna et al.l (|2012l ) hypothesis. 

The non-detection of thermal flux from HAT-P-12b 
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is intriguing since neither the Burrows nor the Fortney 
models can account well for this. Searching the litera- 
ture, we find two other cool low-mass planets o bserved 
with Spitzer at 3.6 and 4 5 /xm - GJ 436b (Butle ret alJ 
[200l iGillon _et al. '2007b'; 'Ste venson et"al.l 12010) and 
WASP-29b (iHardin et al..,2012il . GJ 436b, WASP-29b 
and HAT-P-12b are the three coolest planets observed 
during secondary eclipse with Spitzer. They have equi- 
librium temperatures of ~ 600 K, ^ 980 and ^ 950 K, re- 
spectively, assuming complete redistribution of the stel- 
lar flux to the night side and zero albedo. Assuming 
no flux redistribution and zero albedo, the average equi- 
librium temperatures of their day sides are ~ 700 K, 
~ 1110 K and ^ 1120 K, respectively. For these calcula- 
tions, we assume that the planets are always at a distance 
of 1 semi-major axis away from their host stars, which is 
relevant for GJ 436b, since it has non-zero eccentricity, 
e = 0.16. GJ 436b orbits an M2.5 V star with effectiv e 
temperature, Tes = 3350 ± 300 K (iManess et al.l l2007D. 
and [M/H] = -0.32±0.12 (|Bean et al.ll2006[) . The HAT- 
P-12b and WASP-29b host stars are larger - K4 dwarfs 
with Teff = 4650 ± 60 K, and 4800 ± 150 K, respectively. 
The metallicity of HAT-P-12 is almost identical to that 
of GJ 436: [Fe/H]= -0.29 ± 0.05, but for WASP-29, 
[Fe/H l= 0.11±0.014 ()Hartman et al.|[200l iHellier et al.l 
mm . WASP-29b, like HAT-P-12b has a mass similar 
to Saturn's, while GJ 436b is a hot Neptune. 

WASP-29b and GJ 436b exhibit a measurable eclipse 
depth at 3.6 /im, while none is detected at 4.5 /im. This, 
in combination with Spitzer secondary eclipse measure- 
ments for GJ 436b at 5.8, 8.0, 16 and 24 ^m, prompts 
iStevenson et al.l (|2010[ ) to suggest that the planet has a 
non-inverted atmosphere with large concentrations of CO 
at the expense of CH4. They argue that the small eclipse 
depth at 4.5 /im is a result of a strong absorption feature 
of CO, while the strong 3.6 /im eclipse is caused by lack 
of CH4 absorption. CH4 is expected to begin to domi- 
nate as a carbon bearing molecule below temperatures of 
abou t HOOK (for p ressures o f 1 bar and solar metallic- 
ity) (|Lodders fc Feg lcv 2003 iFortnev et al. 2008), and 
so the suggested abundance of CO would require ther- 
mochemical d isequilibrium in the GJ 436b's atmosphere. 
IHardin et al.l ([21312) suggest that a similar explanation is 
possible for WASP-29b. 

Despite their similar host stars, irradiation levels, plan- 
etary masses and planetary equilibrium temperatures, 
HAT-P-12b exhibits a completely different behavior from 
WASP-29b, since it produces prominent eclipses at nei- 
ther 3.6, nor 4.5 /xm. The reasons for this disparity are 
unclear and additional atmospheric modeling and exten- 
sive observations are needed to investigate the atmo- 
spheres of planets with temperatures < 1000 K, which 
appear to be very dissimilar to those of traditional hot 
Jupiters. 

4.4. Orbital Eccentricity Constraints 

The measured time of secondary eclipse can be used 
to constrain the orbital eccentricity, e, as part of the 
quantity | e cos w | , where uj is the arg ument of periastron. 
For a detailed discussion, see, e.g., iCharbonneau et all 
(|2005D . For HAT-P-3b and HAT-P-4b, we average the 
timed eclipse central phases, weighing them by the in- 
verse of their variance and derive central phases of 
0.50292 ± 0.00076 (HAT-P-3b) and 0.49951 ± 0.00070 



(HAT-P-4b). The individual timings for HAT-P-4b agree 
within 0.1 CT, but the ones for HAT-P-3b disagree at the 
2.8 cr level. We suggest that this discrepancy may be due 
to residual red noise, potentially due to stellar activity 
artifacts in the light curves for this planet that has not 
been taken into account by our data analysis procedures. 

Even for a circular orbit, or for an argument of peri- 
astron of 0/180°, the central phase of the eclipse is not 
expected to be exactly 0.5. The light travel time de- 
lay is 19 s, 22 s and 19 s, corresponding to to expected 
central phases of 0.50008, 0.50008 and 0.50007 for HAT- 
P-3b, HAT-P-4b and HAT-P-12b, respectively. Another 
source of apparent delay could be an off-centered hottest 
point on the planet's face, closer to t he trailing side of the 
planet than to its leading limb (e.g., 'Charbonneau et alj 
2005; Cooper & Showman 2005; Knutson ct al. 2003 
Agol et al.i 2010). Orbital perturbations due to previ- 
ously undetected planets can cause the secondary eclipse 
to occur earlier or later than expected. Our timing pre- 
cision is not sufficient to detect any of these effects. 

Taking the light travel ti me delay into account, ari d 
following the discussion in ICharbonneau et"al] (j2005D . 
we find that for HAT-P-3b, |e cosu;| < 0.0081 within 
3ct. Similarly, within 3cr, |e cosw| < 0.0042 for HAT- 
P-4b. These values are consistent with the value of 
e = adopted by the discovery studies ([Torres et al.l 
I2007t iKovacs et ai.|[2007l ) for all three planets. Since we 
do not detect the eclipses of HAT-P-12b, we cannot con- 
strain its eccentricity. 

5. CONCLUSION 

We measure secondary eclipse depths in the 3.6 and 
4.5 /im bands of Warm Spitzer for three exoplanets, 
HAT-P-3b, HAT-P-4b and HAT-P-12b. We find that 
IIAT-P-3b and IIAT-P-4b have inefficient heat transfer 
from their day to night sides, but the models we compare 
to the data are ambiguous with respect to possible tem- 
perature inversions in their atmospheres. We detect the 
eclipses of HAT-P-12b neither at 4.5 /zm, nor at 3.6 /xm. 
This result is in contrast with Spitzer eclipse measure - 
ments of hot Neptune GJ 436b by Stevenson ct al. (2010) 
and hot Saturn WASP-29b by (jliirdin ct al. 2012), two 
planets in the same temperature regime as HAT-P-12b. 
This is a confirmation that current models need further 
development to explain observations of planets with ef- 
fective temperatures below 1200 K. Additional infra-red 
secondary eclipse observations of "warm Jupiters" are 
urgently needed to constrain these models. Even non- 
detections, as in this study, combined with accurate es- 
timates of the |e cosw| quantity would be invaluable in 
assessing the plausibility of atmospheric models. Cur- 
rently, the Spitzer Space Telescope remains the best avail- 
able observatory to perform these studies. 

In addition, we find that the eccentricities of HAT-P-3b 
and HAT-P-4b are consistent with zero. This is in agree- 
ment with previously published radial velocity measure- 
ments. We compile transit timing measurements from 
the literature and improve the ephemeris of HAT-P-12b. 
We see no obvious transit timing variations that could 
indicate additional objects in the HAT-P-12b system. 

Future eclipse depth observations will benefit from 
the experience that the community has gained from the 
study of exoplanets with the Warm Spitzer mission. The 
pointing oscillation observed for long time series ob- 
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servations with IRAC has been a problem because it 
causes apparent brightness changes of the target stars 
due to intrapixel response variations (the pixel-phase ef- 
fect). This effect has been mitigated, for observations 
after October 17, 2010, by adjusting the operation of 
a heater designed to keep a battery within its operat- 
ing temperature range. Both the amplitude and the pe- 
riod of the pointing wobble have been reduced, resulting 
in smaller brightness changes over shorter periods (from 
60 min to ~ 40 min) , making them more distinguish- 
able from astrophysical phenomena like exoplanet tran- 
sits and eclipses. As a result, the removal of the point- 
ing oscillation effect can be done more efficiently and 
precisely^^. In December 2011, improved pointing using 
the Pointing Calibration and Reference Sensor (PCRS) 
was added to the warm-mission IRAC time series photo- 



metric observations. This option was available and fre- 
quently used during the cryogenic mission and, when it 
is possible to utilize it, can be expected to help reduce 
the pixel-phase effect for staring observations of point 
sources longer than ~ 12hours^ . These improvements 
increase the value of the observatory as one of the best 
currently available tools to study exoplanet atmospheres. 
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TABLE 1 

Adopted Stellar, and Planetary Parameters 



HAT-P-3bd HAT-P-4b° HAT-P-12bg 



M,, (Mq) 


0.917 ±0.030 


1 r,y-| +0.120 

-^■^'-^-0.070 


0.733 ±0.018 


R* (R©) 


0.799 ± 0.039 


1 600+°"'' 
l.DUU_Q 


+0.017 
"JJ^-0.012 


Ks (mag)'" 


9.448 ± 0.025 


9.770 ±0.020 


10.108 ±0.016 


Teff- (K) 


5185 ± 80 


5990f 


4650* 


b imp act 


0.530 ±0.075 




0.211 ± 0.012 


Mp (Mj) 


0.591 ± 0.018 


0.680l°;°r5 


PI r,„ +0.068 


Rp (Rj) 


0.827 ± 0.055 


1 OO7+0.079 

1.0Ci(_Q Q3g 


n 050+0.029 


P (days)'' 


2.8997382 ± 0.0000009 


3.0565254 ± 0.0000012 


3.21305929 ± 0.00000034 


ap (AU)<= 


0.03866 ± 0.00042 


0.04465 ±0.00111 


0.03842 ± 0.00031 



^ Two Micron All Sky Survey (2MASS) Kg magnitude of the star. 

^ The HAT-P-3b and HAT-P-4b orbital periods arc taken Sa da et all ( 120121) : the HAT-P-12b period 
is derived from our updated ephemeridcs shown in Tabic and Figure [2] 

Calculated from the stellar masses and orbital periods assumed in this table. 
^ Values from [Uhan et all ( |2011h , except for the magnitude, Ks, the period, P, and the semimajor 
axis. 

Values from [SouthworthI ||201in . except for the magnitude, Kg, the effective temperature, T^ff , the 
period, P, and the semimajor axis. 
MKnutson et al.| |20lo|). 

S Values from flTartman et al.| j200gt) , except for the magnitude, Kg, the effective temperature, T^ff , 
the period, P, and the semimajor axis. 



TABLE 2 
Photometry Parameters 





HAT 


-P-3b 


HAT 


-P-4b 


HAT- 


P-12b 




3.6 fim 


4.5 ^m 


3.6 fim 


4.5 ^m 


3.6 fj,m 


4.5 ^m 


Tflx (pixels)'' 


2.5 


3.0 


4.0 


3.0 


5.0 


2.5 


Rl (pixels)'' 
b 


3.5 


6.5 


4.0 


6.5 


6.0 


6.5 


0.40 


0.30 


0.85 


0.50 


0.95 


0.50 


c (pixels) 


1.30 


1.60 


0.35 


1.05 


1.05 


1.10 


median rvar 


2.21 


2.33 


2.24 


2.55 


3.10 


2.46 


maximum rvar'* 


2.33 


2.39 


2.47 


2.64 


3.24 


2.58 


minimum rvar'* 


2.16 


2.25 


2.09 


2.41 


2.96 


2.37 



^ Radius adopted for fixed aperture photometry. 

^ Radius adopted for calculating the stellar intensity for Eg nation fl] 
Variable aperture photometry parameters adopted, see Equationm 
Median, maximum and minimum variable aperture radius. 



TABLE 3 
Observation Details 





HAT-P-3b 


HAT-P-4b 


HAT-P-12b 




3.6 fj,m 


4.5 fj,m 


3.6 fim 


4.5 ^m 


3.6 /im 


4.5 fim 


Observation start (UTC) 
Observation end (UTC) 
Orbital phase coverage 
Image count 


17-03-2010,02:45 
17-03-2010,10:28 
0.433 - 0.544 
13,760 


20-03-2010.00:04 
20-03-2010,07:46 
0.428 - 0.539 
13,760 


12-04-2010,02:21 
12-04-2010,09:59 
0.439 - 0.543 
3,871 


02- 09-2010,18:49 

03- 09-2010,02:27 
0.448 - 0.552 
3,871 


16-03-2010,13:58 
16-03-2010,21:35 
0.445 - 0.544 
2,097 


26-03-2010,05:09 
26-03-2010,12:47 
0.443 - 0.542 
2,097 
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TABLE 4 
HAT-P-12B Ephemerides 



N 


Observation Date 


Primary Transit (BJDtt) 


Uncertainty 


Notes 





2007 Mar 28 


2454187.85655 


0.00020 


Hartman et al. (2009) 


9 


2007 Apr 26 


2454216.77265 


0.00014 


Hartman et al. (2009) 


212 


2009 Feb 06 


2454869.02397 


0.00017 


Hartman et al. (2009) 


221 


2009 Mar 07 


2454897.94225 


0.00024 


Hartman et al. (2009) 


238 


2009 May 01 


2454952.56398 


0.00080 


ETD 


242 


2009 May 13 


2454965.41639 


0.00046 


ETD 


242 


2009 May 13 


2454965.41748 


0.00090 


ETD 


248 


2009 Jun 02 


2454984.69368 


0.00060 


ETD 


350 


ZUlU Apr ZO 


z400oiZ.4zD / o 


O.OOOoz 


rj 1 U 


361 


2010 May 31 


2455347.76929 


0.00021 


Sada et al. (2012) 


449 


2011 Mar 10 


2455630.51896 


0.00049 


ETD 


454 


2011 Mar 26 


2455646.58477 


0.00059 


ETD 


454 


2011 Mar 26 


2455646.58486 


0.00040 


ETD 


455 


2011 Mar 29 


2455649.79769 


0.00020 


Lee et al. (2012") 


458 


2011 Apr 07 


2455659.43563 


0.00038 


ETD 


460 


2011 Apr 14 


2455665.86234 


0.00031 


Lee et al. (2012") 


463 


2011 Apr 23 


2455675.49947 


0.00064 


ETD 


464 


2011 Apr 27 


2455678.71445 


0.00045 


Sada et al. (2012") 


469 


2011 May 13 


2455694.78089 


0.00024 


Lee et al. (20121 


472 


2011 May 22 


2455704.42185 


0.00038 


ETD 


These 20 transits yield: 









P = 3.21305929 ± 0.00000034 days 



To = 2454187.85559 ± 0.00011 in BJDtt 

Note: ETD transits were taken from the Exoplanet Transit Database, |http://var2. astro. cz/ETD/^ and compiled by^ee_^t_al^ ^2012). 
For observations after 31 Dec 2008 and before 30 June 2012, TT ^ UTC+66.184s, while for the two 2007 transits 
TT Si UTC+65.184S. 



TABLE 5 
Secondary Eclipse Results 



Eclipse Depth (%) Brightness Temperature (K)'' Eclipse Central Phase BJDtt'' - 2 450 000 O - (min) 



HAT- 


P- 


-3b, 3.6 /xm 


„ -,12+0.015 


iJI J_162 


^r]^'[^+0-00092 

U.dU01d_g QQJJQ 


5272.82936l°:°°264 


01 9+3-8 
zl.Z_4 g 


HAT- 


P- 


-3b, 4.5 /xm 




1268tl5 


50084+° °°^°® 


5275.71660l°°°l°| 


o 9+4.4 
'-'■^-3.0 


HAT- 


P- 


-4b, 3.6 /-tm 


0.i'iZ_(, g-^g 




49945+0 00091 


5298.78653t°;°°2^5 




HAT- 


P- 


-4b, 4.5 /xm 


n 1 99+0.012 

U.1ZZ_Q gj^4 




49960+° °°"° 
u.4yyDU_g QQJQ2 


5442.44368t°;S 




HAT- 


P- 


-12b, 3.6 lira 


< 0.042 


< 970 








HAT- 


P- 


-12b, 4.5^111 


< 0.085 


< 980 









^ The uncertainty of the brightness temperature included here only takes into account the uncertainty of the eclipse depths, but not the uncertainties 
in the stellar properties and the planetary radius. 

^ Time of secondary eclipse central phase, in Baryccntric Julian Date (BJD) based on Terrestrial Time (TT). For these observations, TT 
UTC + 66.184 s, where UTC is the Universal Coordinated Time. 

^ The measured offset from the expected central phase of 0.50008 (HAT-P-3b and HAT-P-4b), with an adjustment for light travel time, (see 
Section |3.1D , in minutes. 




0.44 [US Q.50 0J2 DJ4 0.44 0.^ DM OM 052 0^4 

Phase 

Fig. 1. — The uncorrected time series photometry for HAT-P-3, HAT-P-4 and HAT-P-12 at 3.6 and 4.5 /im during secondary eclipse 
(dots), with the best fit echpse model we obtain (red lines). These models include the astrophysical eclipse, and the instrumental effects: 
a linear ramp and the dependence of the measured intensity on the x and y-position of the stellar image on the detector. All photometric 
points are shown here, including the initial 49min of HAT-P-3 photometry, whic h we later reject. The best fit models only cover data 
that we have adopted for the fits. We discuss the data models in detail in Section 13.21 The HAT-P-12 photometry is shown on a smaller 
scale, since these data have the longest exposure time, and therefore the highest signal to noise per point. We do not detect the secondary 
eclipses of HAT-P-12b, and in these panels their depths are set to zero. 



HAT-P-12b 



c 
U 

6 




100 



200 



300 



400 



500 



N (planet periods) 



Fig. 2. — The difference between the best fit transit times and the observed transit times for HAT-P-12b. The black symbols represent 
transit timings from the discovery paper for this planet (Hartman et al. 2009). Also included are timings from the E xoplanet Transit 
Database ( ETD, http: //var2.astro.cz/ETD/ ) compiled by i Lee et al.i 120121) (red symbols), and the measurements made bv lLee et al.l 120121) 
(blue) and ISada et al.l (2012) (green). The horizontal axis shows the number of periods after the Tq transit. We do not see any obvious 
correlation in the residuals that may suggest a transit timing variation, and hence the presence of an additional companion in the system. 
On the other hand, the number of outliers on this plot may imply that the timing uncertainties in some of these studies are underestimated. 
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3.6 //m 



4.5 



1.004 r HAT-P-3 



1.002 
1.000 
0.998 
0.996 
^ 1.004 

a 

Si 1.002 



0.996 
1.004 
1.002 
1.000 
0.998 
0.996 



HAT-P-3 




Vi#%'i'#Ai' 



0.44 



0.46 



0.48 



0.50 



0.52 



0.54 0.44 

Phase 



0.46 



0.48 0.50 0.52 0.54 



Fig. 3. — The binned time series photometry for HAT-P-3, HAT-P-4 and HAT-P-12 at 3.6 and 4.5 ^m during secondary eclipse after 
correction for instrumental effects (dots). The bin width is 0.0015 in units of orbital phase, corresponding to about 6 min 16 sec (HAP-P-3b), 
6min 36 sec {HAT-P-4b) and 6 min 56 sec (HAT-P-12b). The red lines represent the best fit eclipse model. The best fit coefficients were 
used to correct all data for instrumental effects, but we have only fitted to the data covered by the eclipse model. The HAT-P-12b eclipse 
depths have been set to zero in this plot. 



c 

3 

O 

O 



3.6 fim 



4.5 nm 



^HAT-P-3 ; ; J 


h 1 






' HAT-P-3 1 1 


' 1 T 1 ' 1 ■ ' ' 1 ■ ' ' 1 ■— 
















0.0008 0.0010 0.0012 0.0014 0.0016 0.0018 



0.0008 0.0010 0.0012 0.0014 0.0016 



100 


-HAT-P-12 




-HAT-P-12 




100 


80 






1 1 
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60 


40 










40 


20 










20 
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Eclipse Depth 

Fig. 4. — We show the histograms of eclipse depth values that result from the prayer-bead technique (black) and the MCMC iterations 
(red). The vertical solid lines indicate the best values from the PBMC (black), MCMC (red) and linear regression (blue), which we take 
to be the median of the eclipse depth values of the prayer-bead and MCMC runs. The dashed lines bracket the regions centered on the 
medians that contain 68% of the recorded values for prayer-bead (black) and MCMC (red). Since MCMC and PBMC address different 
noise effects, we adopt as final best values and uncertainties the results from the technique which yields the larger uncertainty region. For 
HAT-P-12, since there is no secondary eclipse detected, the MCMC does not converge, and the prayer-bead is performed with the central 
phase fixed at 0.50007. The blue dashed lines indicate the Itr estimates, based on the linear regression used to determine the eclipse depth 
in the original data sets. We conservatively adopt these as final best values since they are larger than the prayer-bead uncertainties. 
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0.0014 
0.0012 
0.0010 
0.0008 
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0.0004 
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0.0015 
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.ic„=0.1cm' 
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0.0010 
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Fig. 5. — The measured eclipse amplitudes as a function of wavelengths (filled circles) compared to different atmospheric models: a 
blackbo dy planet with a Kurucz m odel for the spectrum of the star (black lines Kurucz 1979), Burro ws planetary atmospheric models (left 
panels, IBurrows et al.ir2007l . 120081 ). and similar atmospheric models by Fortne v et al.i (right panels, 120051 . l2006ai rB. 120081 ) . The downward 
arrows represent 3 rj upper l imits for the HAT-P- 12b eclipse depths. We have assumed effective stellar temperatures for HAT-P-3 from 
IChan et al] 1120111) and from lKnutson et al.l BOIOC ) for HAT-P-4 and HAT-P-12. Left panels: The red lines represent the inverted models, 
with absorption coefficient of the unknown absorber in the upper atmosphere ftabs = O-lj 0.2, and O.lcm^g"^ and heat redistribution 
parameter Pn = 0.1, 0.1, and 0.3 for HAT-P-3b, HAT-P-4b and HAT-P-12b, respectively. The blue lines are models with no inversion 
('^abs = Ocm^g"^) and Pn = 0.3 for all three planets. We over-plot the theoretical eclipse depths resulting from integrating the model 
stellar and planetary fluxes over the IRAC pass-bands (diamonds for the inverted models and triangles for the inverted ones). Right 
panels: In the Fortney model paradigm, where the stratospheric absorbers are assumed to be TiO and VO, HAT-P-3b and HAT-P-12b are 
sufficiently cool, so that the inverted and non-inverted (i.e., containing or not TiO and VO in the upper layers of the atmosphere) models 
are indistinguishable, due to TiO and VO condensation and rain-out. For HAT-P-3b, we show models with / = 0.6 (red line) and / = 0.25 
(blue line); for HAT-P-4b, both models shown have / = 0.5, but the red model has TiO in the upper layers of the atmosphere, while the 
blue one does not. For HAT-P-12b, the red line represents a model with / = 0.6, while the blue line shows a model with / = 0.5. All 
Fortney model atmospheres have solar composition, except the red model for HAT-P-12, which has metallicity 30 times higher than solar. 
For a more detailed explanation on the models see Section 14.11 



